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O . ABSTRACT 

We have analyzed spectra of the Galactic black hole GX 339-4 obtained through simultaneous 76 ksec XMM- 
Newton/EPlC-pn and 10 ksec RXTE observations during a bright phase of its 2002-2003 outburst. An extremely 



oo 



skewed, relativistic Fe Ka emission line and ionized disk reflection spectrum are revealed in these spectra. Self- 
consistent models for the Fe Ka emission line profile and disk reflection spectrum rule-out an inner disk radius 
compatible with a Schwarzschild black hole at more than the 8er level of confidence. The best-fit inner disk radius 
of 2 - 3 r g suggests that GX 339-4 harbors a black hole with a > 0.8 - 0.9 (where r g = GM/c 2 and a = cJ/GM 2 , 
and assuming that reflection in the plunging region is relatively small). This confirms indications for black hole 
spin based on a Chandra spectrum obtained later in the outburst. The emission line and reflection spectrum also 
\ rule-out a standard power-law disk emissivity in GX 339-4; a broken power-law form with enhanced emissivity 

J> ■ inside ~6r s gives improved fits at more than the 8er level of confidence. The extreme red wing of the line and 

£C) ' steep emissivity require a centrally-concentrated source of hard X-rays which can strongly illuminate the inner 

C*~) . disk. Hard X-ray emission from the base of a jet — enhanced by gravitational light bending effects — could create 

the concentrated hard X-ray emission; this process may be related to magnetic connections between the black hole 
and the inner disk. We discuss these results within the context of recent results from analyses of XTE J 1650-500 
and MCG-6-30-15, and models for the inner accretion flow environment around black holes. 

m 

1. INTRODUCTION 2001, Fabian et al. 2003; for BHCs see, e.g., Miller et al. 2002b, 

4^ ■ T ,. f . • j* 1 uv ui i u i u Miller et al. 2004, Miniuti, Fabian, & Miller 2004). 

O . Irradiation of an accretion disk orbiting a black hole by „,._ . . ' . „ • , „ TT ^, 

*— M fujv j a l -r v GX 339-4 is a recurrent, dynamically-constrained BHC 

i a source of hard X-rays can produce a fluorescent Fe Ka ^ . _ IT / „__.. / ... ... 

O. • ,• , • , ■ . u a. t e ,u (Mrh > 5.8 Mpi\ Hynes et al. 2003) in which radio lets 

emission line, which should bear the signatures of the strong ■ , , „ „ , , , , , ,^ „ , 

?H ' i U f* a •<. ,.• i a k-6 /c u- » i inon with v/c > 0.9 have recently been observed (Gallo et al. 

.i— > . Doppler shifts and gravitational redshifts (Fabian et al. 1989; „„„.x ' IT . , . , „; , 

rr, , i /- o ?- u- innn jc .u ui 1 u 1 u 2004). Herein, we report on the time-averaged 76 ksec XMM- 

)Li see also George & Fabian 1991). If the black hole has near- / / % rvMn „ , 

i • / n nno u 7/^**2\ ■ . Newton/EPlC-pn and 10 ksec RXTE spectra of GX 339— 4, ob- 

. . , maximal spin (a ~ 0.998, where a = cJ/GM ), the innermost . ... , . , , , , ^ , • r ^ •. %- • 

>. , , • , u wtc/-^^ a i»i v l i u tamed during a bright phase (near 1 Crab in soft X-rays) of its 

stable circular orbit (ISCO) around the black hole can be as 2002 2003 outburst 
^ \ small as ;•„, = 1.24 >'.. (where r.. = GM/c 2 . nose /•„, = 6 for 



a = 0); this proximity is expected to produce Fe Ka emission 

'c$ • line profiles with strong red wings because of the relative im- 2. OBSERVATION AND DATA REDUCTION 

portance of gravitational red-shifts in comparison to Doppler GX 339 _ 4 wag obseryed whh XMM . Newton for 15 6 kseCi 

shifts (Laor 1991). starting on 29 September 2002 09:06:42 UT (revolution 514). 

In the X-ray spectra of supermassive black holes in AGN The Epic _ pn camera (StrUder et ^ 2m) was operated ^ 



and stellar-mass black holes in Galactic black hole candidates, 



., , , "burst" mode to accommodate the high count rate expected dur- 

skewed Fe Ka line profiles have proved to be extremely impor- • ... , .• rp, „ .. ■ 

„ \ . F , . . . . ,/ .Jl.. ing this observation. The thin optical blocking hlter was used, 

tant diagnostics of the innermost relativistic regime (for AGN, The data were reduced ugi ^ XMM _ Newton suite SAS ver . 

^U e ' B, J TanakaetaL 19 ^ 5; fo l Bl l C t: see ; e ;8" Mlll T eret aL sion 5.4.1, and the guidelines described in the MPE "cookbook" 

2002a; for a review see Reynolds & Nowak 2003) In some (fjee htt p ://waV e. xr ay.mpe.mpg.de/xmm/cookbook). Events 

cases, evidence for black hole spin may be inferred by the line were extracted ^ a strf ^ RAWX (3 j 5 ^ Q 5) yersus RAWy 

shape (for AGN, see, e.g., Iwasawa et al. 1999, Wilms et al. (2 . 5 -178.5) space. Due to the extremely high source flux back- 
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XMM-Newton Observations of GX 339-4 



ground events were not extracted. The events were then fil- 
tered by requiring "FLAG=0" (to reject bad pixels and events 
too close to chip edges) and "PATTERN < 4" (to accept sin- 
gles and doubles), and the spectral channels were grouped by a 
factor of 5 to create a spectrum. The appropriate canned burst 
mode response file was used to fit the spectrum. 

RXTE observed GX 339-4 for 9.6 ksec starting on 29 
September 2002 09:12:11:28. The data were reduced using 
the suite LHEASOFT version 5.2. Standard time filtering (pri- 
marily, filtering-out the SAA) returned net PCA and HEXTE 
exposures of 9.3 ksec and 3.3 ksec, respectively. For this anal- 
ysis, we have only made use of the spectra from PCU-2 (the 
best-calibrated PCU at the time of writing) and HEXTE-A. 
Events from all layers of PCU-2 were combined to make spec- 
tra; "pcabackest" and the bright source background model were 
used to make background spectra. We added 0.75% systematic 
errors to the spectrum from PCU-2 using the tool "grppha". A 
response matrix (combining rmf and arf files) was generated 
using "pcarsp". The HEXTE spectral files were made using the 
standard recipes, and the standard canned responses were used 
to fit the data. The PCU-2 spectrum was fit in the 2.8-25.0 keV 
band (standard for RXTE analysis; see, e.g., Park et al. 2003), 
and the HEXTE spectrum was fit on the 20.0-100.0 keV band 
(the upper bound was fixed by the upper bound over which the 
ionized disk reflection model is valid). These RXTE spectra 
were fit jointly with an overall normalizing constant allowed 
to float between them. An edge was added at 4.78 keV with 
r = 0.1 to correct for an instrumental Xe edge at this energy in 
fits to the PCA spectrum. 

3. ANALYSIS AND RESULTS 

Fits to the XMM-Newton/EPlC-pn spectrum, and joint fits to 
the EPIC-pn and RXTE spectra do not yield formally accept- 
able fits for the models discussed below. A number of sharp, 
narrow calibration uncertainties remain in the EPIC-pn detec- 
tor response. For the most part, these appear as absorption lines 
or edges in the 2-3 keV range (a feature at 2.31 keV appears 
as an emission line), though similar features are present in the 
0.7-2.0 keV band. Most of these features are likely due to 
Au M-shell edges and Si features in the detector; they are re- 
vealed clearly in this observation due to the high signal-to-noise 
achieved. These narrow-band features do not affect measure- 
ments of the Fe Ka emission line profile, but they do affect the 
overall y 2 fit statistic. Due to these complications, within this 
analysis we use the broad-band capacity of RXTE to charac- 
terize the continuum to serve as a guide for fitting the EPIC-pn 
data, and reserve detailed joint fits for future work. All fits were 
made using XSPEC version 11.2.0 (Arnaud 1996). All errors 
in this work are 90% confidence errors. 

The "Laor" model (Laor 1991) describes the line profiles ex- 
pected around a Kerr black hole. The parameters of this model 
are the line energy Eia 0r , the disk emissivity index q (where the 
emissivity is assumed to have a power-law form of J(r) oc r~ q ; 
q = 3is expected for a standard disk), the inner radius of the line 
emission region r,„ in units of r g = GM/c 2 (with a lower limit 
of = 1.235 r g for a = 0.998), the outer line emission radius 
(fixed to the model limit of R ou , = 400 r g ), the disk inclination, 
and the line normalization. 

The "pexriv" reflection model (Magdziarz & Zdziarksi 1995) 
describes reflection from an ionized accretion disk, and does 
not explicitly include Fe Ka line emission. The important pa- 
rameters in this model are the reflection fraction f (f = £l/2w), 
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FIG. 1. — Above: the data/model ratio obtained when the XMM- 
Newton/EPlC-pn spectrum (black) and RXTE/PCA (red) and HEXTE (blue) 
spectra are jointly fit with a model consisting of MCD and power-law compo- 
nents. The 4.0-7.0 keV range was ignored when fitting the data, and the spectra 
were re-binned for clarity. Below: a blow-up of the EPIC-pn spectrum showing 
the extremely skewed line profile. 

the disk temperature, ionization (£ = Lx /nr 2 ), and inclina- 
tion, the power-law folding energy (£/ w), and the index and 
normalization of the power-law flux irradiating the disk. The 
"constant density ionized disk" reflection model (CDID; Ross, 
Fabian, & Young 1999; Ballantyne, Iwasawa, & Fabian 2001) 
is a reflection model which includes Fe Ka line emission, and 
is especially well-suited to high ionization regimes. The rele- 
vant parameters of this model are log(£), /, and the index and 
normalization of the power-law flux irradiating the disk. 

These disk reflection spectra are calculated in the co-rotating 
or fluid frame, and so must be convolved (or, "blurred") with 
the Laor line element describing the Doppler shifts and gravita- 
tional red-shifts expected for an observer in a stationary frame 
at infinity. In blurring the reflected spectra (and, therefore, the 
power-law continuum), we linked the parameters of the emis- 
sion line and blurring function. The disk components were not 
blurred as the multicolor disk blackbody model (MCD; Mitsuda 
et al. 1984) — only an approximation to the Shakura & Sunyaev 
(1973) disk model because it lacks an inner boundary term — 
may be the best-available model for a disk orbiting a spinning 
black hole when torques are present at the inner boundary. 

3.1. Fits to the RXTE Spectra 

In fits to the RXTE spectra, the equivalent neutral hydro- 
gen column density was fixed at Nh = 5.3 x 10 21 cm" 2 (Dickey 
& Lockman 1990) using the "phabs" model. No reasonable 
continuum model gives an acceptable fit to the RXTE spec- 
tra without the addition of a broad Fe Ka emission line. In 
the absence of added line components, the best-fit "canonical" 
MCD plus power-law model gives x 2 /V = 279.8/74, the bulk- 
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motion Comptonization model (Shrader & Titarchuk 1999) 
gives x 2 / l/ = 766.1/74, and an MCD plus CompTT (Titarchuk 
1994) continuum gives X 2 / l/ = 284.9/72 (where v is the num- 
ber of degrees of freedom). Adding Gaussian emission line and 
smeared edge ("smedge") components to the MCD plus power- 
law model as an approximation to a reflection model, a good fit 
is achieved ( x 1 l v = 68.2/68). The parameters measured via 
this model are: kT = 0.84+° f 4 keV, K MCD = 2800 ± 600, T = 
2.5+g-f, K p i = 2.6 ±0.4, E Gauss = 5.8^| keV, FWHM= 3^ keV, 
Kaauss = x 10" 2 , EW= 300^°° e V, E smedge = keV, 
t = 0.2^j |, and W sm = 2+ 2 keV (where K denotes model nor- 
malizations, and the "smedge" energy was constrained to lie in 
the 7.1-9.3 keV range). 

Fits were also made with blurred CDID and pexriv models. 
Blurring the models for the 1.24-400 r g range and assum- 
ing a power-law emissivity of q = 3.0, we found that the data 
are consistent with strong reflection (/ ~ 1) from a highly ion- 
ized inner disk (log(£) = 4.3-4.4). A line equivalent width of 
EW ~ 200 ± 80 eV was meausured using the pexriv model. 

3.2. Fits to the XMM-Newton/EPIC-pn Spectrum 

Due to the hints of an extremely skewed spectrum found with 
RXTE, we began by fitting more physical models to the XMM- 
Newton/EPIC-pn spectrum. Broad line residuals extend down 
to ^3 keV regardless of whether models with additive compo- 
nents or single-component (e.g., bulk motion Comptonization) 
models are used to fit the continuum, which indicates that the 
continuum does not strongly affect the red wing of the line. The 
first model we considered included MCD and power-law con- 
tinuum components, with added Laor line and smeared edge 
components. The power-law index was constrained to lie within 
A(T) < 0.1 of the RXTE-measmed values. With this model, we 
measure the following values: Nh = 5.1 ±0.1 x 10 21 cm" 2 , fcr = 
0.76 ±0.01 keV, K MCD = 2300!^, T = 2.60_ .o 5 , K p , = 2.2^, 
En,,,,- = 6.97_o. 20 keV, q = 5.5*>f, r in = 2. l^j r g , i = 1 1* deg., 
Kiaor = 7.7±^ x 10" 2 , EW^ = 200^ eV, E sm = 7.9!° ' keV, 
t = 0.6^ A V and a smeared edge width of W = 1.0 ± 0.3 keV. 
An inner radius of r,„ = 6 r ? (as per a black hole with a = 0) 
is excluded at more than the 8<r level of confidence, as is a 
standard emissivity index of q = 3. This model gives x 1 l v - 
3456.5/1894; the formally unacceptable fit is due to the cal- 
ibration problems discussed in Section 3.0. A joint fit to the 
XMM-Newton/EPIC-pn and RXTE spectra using these contin- 
uum parameters is shown in Figure 1 . These spectral parame- 
ters, and the timing properties of the source at the time of this 
observation (Homan et al. 2004), indicate that GX 339-4 was 
observed in the "very high" state. 

Using this model, we measure an unabsorbed flux of 2.1 x 
10~ 8 erg cm" 2 s" 1 in the 0.5-10.0 keV band; the power-law 
contributes 35% of this flux. The total line flux is measured 
to be 5.3 x 10~ 10 erg cm" 2 s" 1 , or 8.4 x 10" 2 ph cm" 2 s. It 
should be noted that the best fit without a line component gives 
X 2 /^ = 7552/1902. Clearly, abroad line component is required 
at far more than the 8er level of confidence. 

Models which allow the disk emissivity law to assume a bro- 
ken power-law form, with J(r) oc r"*" within radius r q , and 
J(r) oc r~ q °"' outside of r q , yield statistically improved fits over 
models assuming simple power-law emissivity laws at more 
than the 8er level of confidence. 

Allowing for a broken power-law form for the emissivity and 
a blurred CDID reflection model, we measure: A^/ = 5.5±0.1 x 
10 21 cm" 2 , kT = 0.79 ± 0.01 keV, K MCD = 1860 ± 20, 
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FIG. 2. — The EPIC-pn spectrum of GX 339-4 fit with a relativistically 
blurred reflection model, with data/model ratio shown below. The total spec- 
trum is shown in green, the disk component in blue, and the CDID reflection 
model is shown in red (prior to blurring). The data were re-binned for clarity. 

r = 2.61^-g?, Kcdid = 2.8 ±0.2 x 10" 26 , log(0 = 4.5^, f = 
2.0_o.2, n„ = 2.9 ±0.1, q in = 5.5^, q om = 3.0 +a6 , r, = 7±2r s , 
and i = 1 1^ deg., for x 2 /v = 3887.3/1895. This fit to the EPIC- 
pn spectrum is shown in Figure 2. 

Fits with a blurred pexriv model with a broken power-law 
emissivity (£ = 3.0 x 10 4 , kT^ = 1.0 keV as per Young et al. 
2001, E f id = 200 keV, r q = 6 r g , and q out = 3 were fixed) yielded 
the following fit parameters: Nh = 5.3 ±0.1 x 10 21 cm" 2 , kT = 
0.76 ± 0.01 keV, K MCD = 2200 ± 300, T = 2.61^-g?, K PL . pexriv = 
1.5 ± 0.1, i = \2 + _\ deg., Eiaor = 6.97_o.o7 keV,' r in = 2.1^, 
q in = 5.5^j, Kiaor = 7.2 ± 0.3 x 10" 2 , EW = 210^ eV, and 
/= 1.0!^, giving x 2 /v = 3481.3/1896. 

The measured inner disk color temperature and high disk ion- 
izations are expected in the very high state. Reflection fractions 
greater than unity might be taken as indications of light bending 
effects. All fits show negative residuals above £ ~ 8.0-8.5 keV. 
Reflection models predict a steeper continuum and also account 
for ionized Fe K-shell absorption in this range. The observed 
residuals can be corrected with a phenomenological smeared 
edge component "r" ~ 0.5,but would be better modeled by re- 
finements to reflection models. 

At high mass accretion rates, disks around black holes likely 
extend to the ISCO. Although the Laor model assumes maxi- 
mal black hole spin, it is likely broadly valid over a range of 
high spin parameters. Assuming that r,„ = rjsco (equivalent to 
assuming that the inner disk drags matter in the plunging re- 
gion minimally; see Krolik & Hawley 2002, Reynolds & Begel- 
man 1997, and Young, Ross, & Fabian 1998); equation 2.21 in 
Bardeen, Press, & Teukolsky (1972) may be used to estimate 
the spin parameter of the black hole in GX 339-4. Using this 
procedure, the line and reflection fits detailed above suggest that 
GX 339-4 harbors a black hole with a > 0.8-0.9. 

4. DISCUSSION 

We have discovered an extremely skewed Fe Ka emission 
line in an XMM-Newton/EPIC-pn spectrum of the Galactic 
black hole GX 339-4. Spectral fits with the Laor relativistic 
disk line model and a phenomenological continuum, and with 
the Laor model and self-consistent blurred disk reflection mod- 
els, provide tight constraints on the nature of the black hole and 
inner accretion flow in GX 339-4. Our fits indicate that the in- 
ner disk likely extends to r in = 2-3 r g , translating to a black 
hole spin parameter of a > 0.8-0.9. Moreover, an enhanced 
inner disk emissivity index of q = 4.8-6.0 is required within 
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5-9 r g . Fits with the same models which fix either r,„ = 6.0 r g 
(as per an a = Schwarzschild hole) or q = 3.0 (as per a standard 
disk), are more than 8(7 worse than the best-fit values summa- 
rized here and detailed in the previous section. 

X-ray spectroscopy of XTE J 1650-500 (Miller et al. 2002b, 
Miniutti, Fabian, & Miller 2003) and MCG-6-30-15 (Wilms 
et al. 2001, Fabian et al. 2002) has revealed remarkably sim- 
ilar line profiles and inner disk emissivity properties; both 
XTE J 1650-500 and MCG-6-30-15 may have spin parameters 
which are similar to GX 339-4. Comparisons of the timing 
phenomena observed in stellar-mass Galactic black holes and 
Seyfert galaxies have shown that noise properties may simply 
scale with mass (Uttley & McHardy 2001). It has also recently 
been shown that Galactic black holes may have AGN-like warm 
absorbers (Miller et al. 2004). Our results extend these findings 
in that they suggest that — at least in certain states — details 
like enhanced inner disk emissivity are likely similar. 

The models we have considered not only constrain the black 
hole spin parameter, but also the nature of the inner accretion 
flow geometry. The inner disk must be strongly illuminated 
by a centrally-concentrated source of hard X-rays. Miniutti & 
Fabian (2003) have calculated the effects of light bending on 
Fe Ka line variability. This model can explain the complex 
time variability seen in MCG-6-30-15, and appears to describe 
the behavior of the line strength in XTE J 1650-500 remarkably 
well (Miniutti, Fabian, & Miller 2004; Rossi et al. 2003). The 
light bending model assumes a power-law source of hard X-ray 
emission above the black hole which moves "vertically" along 
the black hole and inner disk angular momentum axis. While 



models for jets in Galactic black holes require that X-ray emis- 
sion be focused away from the disk (see, e.g., Markoff, Falcke, 
& Fender 2001), at the base of the jet the flow is less relativis- 
tic and the beaming therefore less extreme. It is possible that 
synchrotron self-Compton emission from the base of a jet may 
be focused onto the inner disk by light bending (this would rep- 
resent an extension to the jet reflection considered by Markoff 
& Nowak 2004). This does not rule-out additional hard X-ray 
emission from a corona. 

Heightened inner disk emissivity of the kind seen in 
GX 339-4 can plausibly be explained by the dissipation of the 
black hole's rotational energy via magnetic connections to the 
inner disk (Blandford & Znajek 1977; see also Gammie 1999 
and Li 2003). The emissivity we have measured is above that 
predicted by magnetic connections to matter in the plunging re- 
gion (q = 7 /2, Agol & Krolik 2000), though it is possible this 
process could also be at work. It is possible that such processes 
and the putative illumination of the inner disk region by the 
base of a jet might be intimately related; indeed, the extraction 
of black hole spin energy is often invoked as a possible means 
of powering relativistic jets in AGN. 
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